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Abstract. Let [c, d] be an interval on the real line and (i be a measure of 
the form d/j, — fj,duji c<< j\ with fi = hh, where h(t) = (t — c) ac (d — t) ad , 
a c ,ad £ [0,1/2), h is a Dim-continuous non-vanishing function on [c, d] 
with an argument of bounded variation, and W[ c d j is the normalized arc- 
sine distribution on [c, d] . Further, let p and q be two polynomials such that 
deg(p) < deg(g) and [c, d] Pi z(q) = 0, where z(q) is the set of the zeros of q. 
We show that AAK-type meromorphic as well as diagonal multipoint Pade 
approximants to 



converge locally uniformly to f in Q\ n D and £>(, respectively, where Df is 
the domain of analyticity of f and D is the unit disk. In the case of Pade 
approximants we need to assume that the interpolation scheme is "nearly" 
conjugate-symmetric. A noteworthy feature of this case is that we also allow 
the density /i to vanish on (c, d), even though in a strictly controlled manner. 
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1. Introduction 

Let f be a function of the form 



where [c, d] = supp(/z) is the support of a complex Borel measure /i, the polyno- 
mials p and q are coprime, deg(p) < deg(q) —: m, and z(q) is the set of zeros of q. 
Let iO[c,d] be the equilibrium distribution for [c,d], which is simply the normalized 
arcsine distribution. In this paper, we assume that \i is absolutely continuous with 
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respect to W[ C)( q and fi, its Radon-Nikodym derivative (dfj, = p,du)[ c ^), is such that 

fi = hh, (1.2a) 
ij. = hhh Xl (1.2b) 

where ft, is a non-vanishing Dini-continuous function with argument of bounded 
variation on [c,d], H(t) :— \t — c\ ac \t — d\ ad , a c ,ad 6 [0, 1/2), x C (c,d) is a finite 
set of distinct points, and h x (t) :— Y\ xex \t — x| 2ax , a x € (0,1/2). Under such 
assumptions on f, we show locally uniform convergence of L p (T)-best meromorphic 
(in this case we assume that [c, d] C (—1, 1)) and certain diagonal multipoint Pade 
approximants to f in 

3), :=C\(supp( A1 )Uz(q)), (1.3) 

the domain of analyticity of f, where C is the extended complex plane. It is known 
[25l |5] that the denominators of both types of approximants satisfy non-Hermitian 
orthogonality relations with respect to /x that assume a similar form. This leads to 
similar integral representations for the error of approximation, which is the reason 
why we treat them simultaneously. 

Generally speaking, meromorphic approximants (MAs) are functions mero- 
morphic in the unit disk that provide an optimal approximation to f on the unit cir- 
cle in the L p -norm when the number of poles is fixed. When considering them, it is 
customary to assume that supp(/x)Uz(<7) is contained in the unit disk, D. The study 
of MAs originated from the work of V.M. Adamyan, D.Z. Arov, and M.G. Krein 
PP, where the case p = oo was considered. Nowadays such approximants are often 
called AAK approximants. The L p -extensions of the AAK theory were obtained 
independently by L. Baratchart and F. Seyfert [5] and V.A. Prokhorov [21J . Mero- 
morphic approximation problems have natural extension to Jordan domains with 
rectifiable boundary when the approximated function f is meromorphic outside of 
a closed hyperbolic arc of this domain [3]. However, we shall not consider such a 
generalization here. 

The AAK theory itself as well as its generalizations is based on the intimate 
relation between best (locally best) MAs and Hankel operator whose symbol is the 
approximated function [UG3I5T]. The study of the asymptotic behavior of MAs is, 
in fact, equivalent to the study of the asymptotic behavior of the singular vectors 
and singular numbers of the underlying Hankel operator (see Section [3]). Hence, 
the present work (more specifically, Theorems [I] and [2]) can be considered as an 
asymptotic analysis of the singular vectors of Hankel operators with symbols of 
type ([Ll|-( |1.2a| . 



Let us briefly account for the existing results on convergence of MAs to 
functions of type Uniform convergence was obtained in [3] for the case p = 2 

(in this case meromorphic approximants reduce to rational functions) whenever 
/i is a positive measure and the rational summand is not present, i.e. q = 1 and 
necessarily p = 0, (such functions f are called Markov functions). The general case 
p G [1, oo] was addressed in [1], where again only Markov functions were considered 
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and uniform convergence was shown under the assumption that u belongs to the 
Szego class, i.e. log(dn(t) / dt) is integrable on [c, d]. The case of complex measures 
and non-trivial rational part was taken up in [5], where convergence in capacity 
in D \ supp(/x) was obtained while supp(/i) was assumed to be a regular set with 
respect to the Dirichlet problem and /i had to be sufficiently "thick" on its support 
and have an argument of bounded variation. 

On the other hand, diagonal multipoint Padc approximants (PAs) are ratio- 
nal functions of type (n, n) that interpolate f in a system of In not necessarily 
distinct nor finite points (interpolation scheme) lying in S)j with one additional 
interpolation condition at infinity. Unlike the meromorphic case, it is pointless to 
assume that supp(/x) and z(q) lie in D. It is customary to call PA classical if all 
the interpolation points lie at infinity. Such approximants were initially studied by 
A. A. Markov [17] using the language of continued fractions. Later, A. A. Gonchar 
[13] considered classical PAs to functions of type ( |1.1[ ) with nontrivial rational 
part and positive u. Locally uniform convergence to f in D j was obtained un- 
der the condition that fi belongs to the Szego class. Continuing this work, E.A. 
Rakhmanov has shown [35] that the restriction on /j to be in the Szego class 
cannot be relaxed in general, but if all the coefficients of R are real, uniform con- 
vergence holds for any positive measure. In the recent paper [14] A. A. Gonchar 
and S.P. Suetin proved that uniform convergence of classical PAs still holds if u is 
a complex measure of the form dfj, = hduj^ c d ^ , where h is a non- vanishing analytic 
function in some neighborhood of [c, d\. Recently, using the operator-theoretic ap- 
proach, M.S. Derevyagin and V.A. Derkach [TT] showed that there always exists 
a subsequence of diagonal PAs that converges locally uniformly to f whenever the 
latter is such that /i is a positive measure and p/q is real-valued on supp(/x) but 
can have poles there. Finally, we mention a weaker result that holds for a larger 
class of complex measures. It was shown in [5] Thm. 2.3] that multipoint Pade ap- 
proximants corresponding to "nearly" conjugate-symmetric interpolation schemes 
converge in capacity in C \ supp(/i) whenever supp(/z) is a regular set with respect 
to the Dirichlet problem and u is sufficiently "thick" on its support and has an 
argument of bounded variation. 

The main results of this paper are presented in Section [3] Theorems [T]^4j and 
Section [4] Theorems [5] and [6] The conditions imposed on the measure \i in these 
theorems come from Theorem [7] The latter is, in fact, a consequence of Theorems 
2 and 3 in [7J. In particular, if Theorem |jj is established under other assumptions 
on /i, this would yield Theorems [T}j6] for this new class of measures. For instance, 
all the main results of the present work would hold whenever fj, is of the form 

du(t) = h(t)(t - c) a "(d - t) ad dt, 

where h is an m-times continuously differcntiable non- vanishing function on [c, d] 
with m-th derivative being ^-Holder continuous and a c ,ctd <G (— l,oo) PI (—to — 
<t,to + |S]. 
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2. Preliminaries and Notation 

To smoothen the exposition of the material in main Sections [3] and [4] we gather 
below some necessary prerequisites and notation. 

Let T s := {z : \z\ = s}, Tf := T s (l{z : ±Im(z) > 0}, and D s := {z : \z\ < s}, 
s > 0, be the circle, the semicircles, and the open disk centered at the origin of 
radius s. For simplicity, we drop the lower index 1 for the unit circle (semicircles) 
and the unit disk. 

Denote by H p , p G [l,oo], the Hardy spaces of the unit disk consisting of 
holomorphic functions / in D such that 

11/11?:= sup ~ / \f(s0\ p m <oo if pe[l,oo), 

o< s <i zir J T ^2.1) 

|| /Hoc := sup |/(^)| < oo if p = oo. 

It is known [T^l Thm. 1.5.3] that a function in H p is uniquely determined by its 
trace (non-tangential limit) on the unit circle and that the L p -norm of this trace 
is equal to the iJ p -norm of the function, where LP is the space of p-summable 
functions on T. This way H p can be regarded as a closed subspace of LP . 

In the same vein, we define Hq, p € [1 , o o] , c onsisting of holomorphic func- 



tions in C\D that vanish at infinity and satisfy (2.1 ) with 1 < s < oo and z € 



respectively. In particular, we have that L 2, = LF 1 © Hq. Thus, we may define or- 
thogonal projections V+ : L 2 — > H 2 (analytic) and V- : L 2 — * Hq (antianalytic). 
It is easy to see that 

MO d£ _ f V+(h)(z), zeD, h( _ L2 



T £-z27ni [ -V_(h)(z), zeC\D 

Recall also the well-known fact [T2l Cor. II. 5. 8] that any nonzero function in H p 
can be uniquely factored as h = jw, where 

^)=exp{^||^log|/i(0IK|}, zeD, 

belongs to H p and is called the outer factor of h, while j has modulus 1 a.e. on T 
and is called the inner factor of h. The latter may be further decompose as j = bs, 
where 6 is a Blaschke product, i.e. a function of the form 



\zA 1 — ZjZ 



that has the same zeroing as h, while s is the singular inner factor. For simplicity, 
we often say that a function is outer (resp. inner) if it is equal to its outer (resp. 
inner) factor. 

Continuing with the notation, for any point-set K and any function / G H p , 
we denote by K* and f a their reflections across T, i.e., K* :={z:l/zC K} and 
f a (z) :— z^ 1 f(l/z). Clearly then / £ Hq and the map is idempotent. Further, 
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for an interval [c, d] we set 

«:=4/(d-c), t»(z) := y/{z - c)(z - d), and D :=C\[c,d], (2.2) 

where such a branch of ro is chosen that in is holomorphic in D \{ 00} and w(z)/z — > 
1 as z — * 00. Then 

5(«) := l/(l/tt)(z)) CT = z 2 ro CT (z) = ^(1 - cz)(l - dz), 

is holomorphic in D* \ {00} and fo(0) > 0. Moreover, the function 

2z- (d+c) - 2vo(z) _ 

i>{z):= K — -1 K -L, zeT), 2.3 

a — c 

is the conformal map of J) onto D such that ip(oo) = and ip'(oo) = k > 0. It is 
also easy to see that ifi has well-defined unrestricted boundary values from both 
side of [c, d] (we assume that the positive side of [c, d] lies on the left when the 
interval is traversed in the positive direction, i.e. from c to d). Moreover, it holds 
that 

= 1 on [c,d]. (2.4) 
Let now h be a Dini-continuous non-vanishing complex-valued function on 
[c, d] . Recall that Dini-continuity means 

S^ 1 max \h(ti) — h(t2)\dS < 00. 

[0,d-c] *a|<« 

It can be easily checked (cf. [3 Sec. 3.3]) that the geometric mean of ft., i.e. 

G h := exp |y log /i(t)do; [C;d] (t) J , (2.5) 

is independent of the actual choice of the branch of the logarithm and is non-zero. 
Moreover, the Szego function of h, i.e. 

e*pj^ j l °^Y du ^-\ J logWW*)}' ( 2 - 6 ) 

z G ©, also does not depend on the choice of the branch (as long as the same 
branch is taken in both integrals) and is a non- vanishing holomorphic function in 
D that has continuous boundary values from each side of [c, d] and satisfies 

h=G h S+S~ on [c,d] and S h (oo) = 1. (2.7) 

The continuity of the traces of Sh is ensured by the Dini-continuity of h, essentially 
because Dini-continuous functions have continuous conjugates [12, Thm. III. 1.3]. 
In fact, more can be said. Let a be a non- vanishing holomorphic function in D 
that has continuous traces on each side of [c, d] and a(oo) = 1. Suppose also that 
a + a~ = c 2 for some constant c. Then the functions dioip ;= c/a and a e o (X/ip) := 
a/c are holomorphic in D and C \ D, respectively, have continuous traces on T, 
Oi(0) = c and a e (oo) = 1/c. Moreover, it can be readily verified that the traces of 
at and a e coincide. Thus, ai and a e are analytic continuations of each other, from 
which we deduce by Liouvillc's theorem that c = 1 and <Zj = a e = 1. This simple 
observation implies the following. Let /i be a Dini-continuous function on [c, d] 
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and G be some constant. If S is a non- vanishing holomorphic function in D that 
assumes value 1 at inifnity, has continuous traces, and is such that GS~ S + = h 



then necessarily G = Gh and S = Sh- It is also true that (2.6) is well-defined 



whenever h is a non-negative integrable function with integrable logarithm; like, 



for example, h and defined after (1.2 1 



We also emphasize that the Szego function of a polynomial can be computed 
in a rather explicit manner as we will now see. Let v be a polynomial with zeros 
in 2), deg(w) < k. Set 

e tU vi-^#(e); 

where z(v) is the set of zeros of v and 771(e) is the multiplicity of e G z(v). Then 
rk(v; •) is a holomorphic function in £) with a zero of multiplicity m(e) at each 
e G z(v) and a (possible) zero of multiplicity k — deg(7j) at infinity. Moreover, it 
has unrestricted continuous boundary values from both sides of [c, d] such that 

r+(v;-)r-(v;-) = l (2.9) 



by (2.4). Then since S v is the unique function of Szego type such that S^S V is 



equal to a constant multiple of v, it holds that 

In some cases it will be important to consider the ratio of the boundary values 
of Szego functions rather then their product. Hence, we introduce 

Q^(t):=S£(t)/S£(t), t€[c,d]. (2.11) 

When h is non- vanishing Dini-continuous function, are continuous on [c, d] and 
assume the value 1 at the endpoints. Finalizing the discussion on Szego functions, 
let us state two of their properties that we shall use implicitly on several occasions 
and the reader will have no difficulty to verify. The first one is the multiplicativity 
property, i.e. Sh 1 h 2 — S^S^, and the second one is the convergence property 
which says that Sh n — [1 + o(l)]S'/ l uniformly in C, i.e. including the boundary 
values, whenever h n = [1 + o(l)]h uniformly on [c,d]. For more information on 
Szego functions of complex h, the reader may consult Sec. 3.3]. 

Next, we denote by A SljS2 := {z : s± < \z\ < S2}, < Si < S2, and 
A s := A s l / S , s < 1, the annuli centered at the origin and by tp the conformal 
map from Q n D* onto A p , ip(l) = 1. Recall that annuli A s are not conformally 
equivalent for different s and therefore p = p([c,d\) is uniquely determined by 
[c, d] . From the potential-theoretic point of view p can be expressed as 



CX ^- cap( M ],T) .^ (2 ' 12) 
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where cap([et, b],T) is the capacity of the condenser ([c, d],T). The map ip is given 
by [241 Thm. VIII.6.1] 



^(z) = exp|r 2 ^ 



dt 
(»»)(*) 

with integration taken along any path in 2) n D* , where 



7 := 



fix 



Mo,i] v/(l-a; 2 )((l-cd) 2 -(d-c) 2 a: 2 )' 



(2.13) 



(2.14) 



Moreover, it holds that tp(z) = tp(z) and <p(l/z) = l/tp(z), z e S n ©*. Thus, 
y(T) = T and </?(B \ [c, d]) = A p x . Further, it is not hard to check that ip ex- 
tends continuously on each side of [c,d] (resp. [c, d]*) and (/j ± ([c, d]) = T p (resp. 
ip ± ( [c, d] * ) = Finally, the Green equilibrium distribution (which is a proba- 

bility measure on [c, d] [24j Sec. II. 5]) for the condenser D \ [c, d], as well as for the 
condenser ID fl Q*, is given by 



dU([c,d],T)(t) 



T 2 dt 



WHt)\ 

Tip 



(2.15) 



where the normalization follows from (2.141 and the second equality holds by 
differentiating (2.131 and taking boundary values. 



3. Meromorphic Approximation 

The meromorphic approximants (MAs) that we deal with arc defined as follows. 
For p£ [1, oo] and n G N, the class of meromorphic functions of degree n in L' p is 

HI := HVB-\ (3.1) 

where B n is the set of Blaschke products of degree at most n (with at most n 
zeros). By the celebrated theorem of Adamyan, Arov, and Krein [T] (see also [501 
Ch. 4]) and its generalizations [ST] it is known that for any fixed n £ N and 
p € [1, oo] and given /eL p there exists a meromorphic function g n such that 

||/-S„|| P = inf 11/ -<?||p. (3.2) 

Moreover, g n is unique when p € [l,oo), but in the case p = oo it is necessary to 
assume / 6 H°° + C(T) to ensure uniqueness of g n , where C(T) is the space of 
continuous functions on the unit circle. Obviously, when supp(/Li) C D and q has 



no zeros on T the function f defined in (1.1 1 complies with these requirements for 



any p 6 [1, oo]. When p < 2, no functional representation for the error is known to 
satisfy orthogonality relations [5]. This is the reason why in what follows we shall 
restrict to the case p € [2, oo]. 
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Due to similar functional decomposition and their appearances in the compu- 
tation^] we consider not only best MAs but more generally critical point of mero- 
morphic approximation problem (3.2). Although their definition is rather technical 
(see below), critical points are just those g n — h n /b n G (see (3.1 1) for which 
the derivative of ||/ — g n \\ p with respect to b n G B n and h n G H p does van- 
ish [5J. By definition, a function g n is a critical point of order n in meromorphic 
approximation problem (3.2) if and only if it assumes the form 

9n = / _«zM = ?MM, (3. 3) 

where Tlf is Hankel operator with a symbol /, i.e. 

Hf.HP'^Hl H f (h):=P-(fh), l/p+l/p' = l/2, 

and v n G H p is of unit norm (a Blaschke product if p = 2), its inner factor is a 
Blaschke product of exact degree n, and is such that 

H}H f (v n ) = a 2 n V+ (\v n \P'~ 2 v n ) if p > 2, 
H}H f (v n ) = V+ (\Hf(v n )\ 2 v n ) if p = 2, 

with TL*j being the adjoint operator. A function v n is called a singular vector 
associated to a critical point g n and 

°n ■= 11/ - 9n\\p = \\Hf(Vn)\\p, pG[2,Oo], (3.4) 

is called the critical value associated to g n . In the case when g n is a best MA to / 
it also holds that a n is the n-th singular number of Tlf, i.e. 

oYi = cr n (7if) := inf — IT 1 1 : T : H p — » Hq linear operator of rank < nX ; 

when p — 2 it is assumed in addition that T is weak* continuous. Hereafter, we 
use the following notation for the inner-outer decomposition of singular vectors: 



v n = b n W n , b n = q n /q n , w n (0) > 0, q n (z) = z n q n (l/z), (3.5) 

where w n is an outer factor, q n is a monic polynomial of exact degree n, and is 
the reciprocal polynomial of To uniformize the notation, we simply set w n = 1 
when p = 2. 

A critical point of order n may have less than n poles, even though we insisted 
in the definition that v n has exactly n zeros. Cancellation may occur due to zeros 
ofV+(fv n ). When this is not the case, we shall call g n an irreducible critical point. 
It is worth mentioning that when p G [2, oo) a best MA is not necessarily unique, 
but has exactly n poles. Thus, all best MAs are irreducible critical points. To the 
contrary, if p — oo, best MA is unique and is the only critical point of order n, but 
may have less then n poles. However, there always exists a subsequence of natural 
numbers for which best AAK approximants are irreducible. Since the behavior of 
the poles of MAs is entirely characterized by this subsequence, hereafter we say 



x It is most likely that a numerical search ends up with stable critical points, that is locally best 
MAs, rather than just best MAs. 
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"a sequence of irreducible critical points" to mean if p = oo that we pass to a 
subsequence if needed. 

Now, we are ready to state the first theorem of this section. 



Theorem 1. Let f be given by (1.1) and (1.2a I. Further, let {g n } be a sequence 
of irreducible critical points of the meromorphic approximation problem to f, p € 
(2,oo]. Then the outer factors w n in (3.51 are such that 

<72 = r + o(i) + ^ ro _^ 

ro q 



where o(l) holds locally uniformly in £>*, T was defined in (2.141, and the poly- 
nomials l n , deg(Z„) < m, converge to zero and are coprime with q. 

This theorem is a strengthening of Lemma 3.4 in [3] that asserts, under much 
milder assumptions on fi, that is a normal family in Sj" and any limit 

point of {w n } in D is zero free. 

For simplicity, set w := (T/ro) 2 ^ for each p > 2. Then Theorem [l] yields 
that w n — > w, uniformly in some neighborhood of D (it follows from the proof of 
Theorem[l]and can be seen from asymptotic formula ( |3.6[ ) that the outer factors w n 
can be extended to holomorphic functions in any simply connected neighborhood 
of D contained in 55* ). Now, we are ready to describe the asymptotic behavior of 
irreducible critical points. 

Theorem 2. Let f and {g„} be as in Theorem^ Then the Blaschke products b„ in 
(3.5 1 are such that 

b n = [1 + o{l)]bip n ~ m D~ 1 locally uniformly in D f nSj, (3.7) 

where b := q/q, m — deg(q), {D n } is a normal family of non-vanishing functions 
in ID fl D* , such that |-D„ | are uniformly bounded above and away from zero on 
[c,d] and [c,d]*. Moreover, the following error estimates take place 

°n - ||f-fl„|| P - [2A/T- 1 + o(l)] P 2 («"™), (3.8) 



where a n is the critical value associated to g n via (3.4|, and 

wj^m*-^ (3 . 9) 
wro \<p J b 1 

uniformly on compact subsets of Df D D, where 

M :=expH log|(6 2 u./ t )(t)|dw (M]!T) (i)|, (3.10) 

is the geometric mean of \b 2 wfi\ with respect to the condenser 2) (~l D* . 

It is worth mentioning that the functions D n are, in fact, Szego functions 
for the condenser J) D D* that first were introduced in [6l Def. 2.38] for the case 
of a positive measure /z. In such a situation the Szego function for a condenser 
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has an integral representation that is no longer valid for complex measures. More- 
over, the normalization in the complex case is more intricate (see Proposition 



111. Nevertheless, it still holds that the functions D n have zero winding number 



on any curve separating [c,d] from [c, d]*, D n (z)D n (l/ z) = 1, z G ID D ID*, and 
M\D+D~\ = \b 2 wfi\ on [c,d]. 

We remind the reader that the case p — 2 has a couple of special traits. First, 
best MA g n specializes to a rational function. Indeed, g n can be written as a sum 
K +_p n -i/q n , where h n G H 2 , deg(p„_i) < deg(q n ) = n. As L 2 = H 2 © H§ and 
f G Hq, we have that 



lf-<?lll = IHl2 + llf-P»-l/<?nll 



Hence, to achieve the minimum of the left-hand side of the equality above, one 
necessarily should take h = 0. This is the reason why we referred on some occa- 
sions to the meromorphic approximation problem with p = 2 as to the rational 



approximation problem. Second, the outer factors w n in (3.5 1 are not present, or 



better assumed to be identically 1. The latter allows us to consider a slightly larger 



class of measures, namely those given by (1.2b I 



Theorem 3. Let f be given by (1.1) and ( 1.2b| ), where sin(a x 7r) G (O^x) 



:= min { lQ^(,)l}exp ^ j , a e x> 

Vh is the total variation of the argument of h on [c,d], sq := maxxlV'l) and 
si := max.j\tp'\. Let further {g n } be a sequence of irreducible critical points of 
the meromorphic approximation problem with p = 2 to f. Then (3.6), (3.7), and 



(3.8 1 hold with w = 1. 



It follows from (3.7) that each b n has exactly m zeros approaching the zeros 



of b. In fact, it is possible to say more. 

Theorem 4. For each r] G z(q) and all n large enough, there exists an arrangement 
of T]i,ni ■ ■ ■ r lm(-q).n, the zeros of b n approaching r\, such that 

where the sequences {maxj. 1^ „|} an d {niax^ (l/A^ |} are bounded above. 

This theorem essentially says that each pole r\ of f attracts exactly m{r\) poles 
of g n , the latter converge geometrically fast and are asymptotically distributed as 
the roots of unity of order m(r]). The proof of this theorem is an adaptation of the 
technique developed in [14 for classical Pade approximants to Cauchy transforms 
of analytic densities. As one can see from the next section, similar results hold not 
only for classical but more generally for multipoint Pade approximants to Cauchy 
transforms of less regular measures. 
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4. Multipoint Pade Approximation 

Let f be given by ( |1.1[ ). Classically, diagonal (multipoint) Pade approximants to 
f are rational functions of type (n, n) that interpolate f at a prescribed system of 



In + 1 points. However, when the approximated function is of the form (1.1 1, it is 
customary to place at least one interpolation condition at infinity. More precisely, 
let <§ = {E n } be a sequence of sets each consisting of 2n not necessarily distinct 
nor finite points in (interpolation scheme), and let v n be the monic polynomial 
with zeros at the finite points of E n . 



Definition (Pade Approximants). Given f of type {1.1) and an interpolation scheme 
$ , the n-th diagonal Pade approximant to f associated with S is the unique rational 
function Yl n = p n /q n satisfying: 

• degp„ < n, deg q n < n, and q n ^ 0; 

• (q n (z)f(z) —p n (z)) /v n (z) is analytic in £> f ; 

• (q n (z)j(z) -p n (z)) /v n (z) = O (l/z n+1 ) asz^oo. 

A Pade approximant always exists since the conditions for p n and q n amount 
to solving a system of In + 1 homogeneous linear equations with 2n + 2 unknown 
coefficients, no solution of which can be such that q n = (we may thus assume 
that q n is monic); note the required interpolation at infinity is entailed by the last 
condition and therefore U n is, in fact, of type (n — l,n). 

By the very definition, the behavior of Tl n depends on the choice of the inter- 
polation scheme. We define the support of § = {E n } as supp(tf) := n ne tqUk>nEk- 
Hereafter, the counting measure of a finite set is a probability measure that has 
equal mass at each point counting multiplicities and the weak* topology is un- 
derstood with respect to the duality between complex measures and continuous 
functions with compact support in C. 

Definition (Admissibility). An interpolation scheme S is called admissible if 

• there exist rearrangements A n of E n such that the sums ~^2 eeE — 
V>(A„(e))| are uniformly bounded when n — > oo; 

• supp((?) C Sf and the probability counting measures of points in E n converge 
weak* to some Borel measure with finite Green energ^\ relative to D . 

Then the following result holds. 

Theorem 5. Let {n„} be a sequence of diagonal Pade approximants associated with 



an admissible interpolation scheme S = {E n } to f given by (1.1 1 and ( |1.2b I with 
a x ii e (O.arcsinTa;), := li: 

for any x G x, where r n := r2 n {v n ', ■)■ Then 



a x n £ (0,arcsinT K ), := liminf min {|(r„Q/ l ) ± (j:)|} , (4.1) 

n— >oo ± 



(f-n„)tt> = [2G A + o(l)] Sj 5 (4.2) 



2 For information on the notions of potential theory we refer the reader to the monographs .23 , 24 . 
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locally uniformly in Dj, where r :— r m (q; ■). When (p/q) = it is not necessary to 
assume boundedness of the variation of the argument of h. 

We would like to point out that is, in fact, continuous function of x on 
[c,d] such that T c = = 1 and min^g^ ^ T x > 0. The latter is true since the 
functions are non- vanishing and continuous on [c, d]. Moreover, it will be shown 
in the proof of Theorem [7] that the admissibility of S implies uniform boundedness 



of |r„ | and hence their uniform boundedness away from zero by (2.9 1. It is also 
easy to check that when the sets E n are conjugate-symmetric and h is a positive 
function, it holds that = 1. 

Concerning the behavior of II„ near polar singularities of f, i.e. near z(g), the 
following theorem asserts the same "roots of unity" behavior as in Theorem [3] and 
is a generalization of [HI Thm. 3] for the case of multipoint Pade approximants 
and less regular measures. 

Theorem 6. Under the conditions of Theorem^let q n be the denominators ofH n . 
Then 

q n = u n ^ m q n<m and q n>m = (1 + o(l))q, (4.3) 

where deg(u„ have no zeros 

on any closed set in £> for all n large enough, and o(l) holds locally uniformly in 
C\z(q). Moreover, for each rj 6 z(g) with multiplicity m(rf) and all n large enough 
there exists an arrangement of rj\ >n , . . . , Tj m (yft n € z(g„ jrra ) such that 

Vk,n = V + Al n [r n ( V )] 1/mM exp i^X , k=l,..., mfa), (4.4) 

where the sequence {maxt 1-^2 n |} * s bounded above. 



5. Non-Hermitian Orthogonal Polynomials 

In this section we describe the asymptotic behavior of non-Hermitian orthogonal 
polynomials with varying weights on [c,d]. In what follows, we assume that {i> n } 
is a sequence of complex measures on [c, d] such that 

dv n = u n duj[ c . d ], v n = hh n hh x ./v n+m , m € Z+, 

where h is a non- vanishing Dini-continuous function on [c, d], {h n } is a normal 
family of non- vanishing functions in some neighborhood of [c, d] none of which limit 
points can vanish in this neighborhood, v n , deg(w„) < 2n, are monic polynomials 
with zeros at finite points of an admissible interpolation scheme, h and h\ are as 
in the introduction with 

a x Tr G ((^arcsinTz), := liminf min {|(r n+m (3/ l ; i?l ) ± (x)|} , (5.1) 

n— >oo ± 

for each x G x, where r*fc := r2k{vk] •)• Observe that we do not require h to have 
argument of bounded variation. Then the following theorem holds. 
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Theorem 7. Let {v n } be as described and {u n } be a sequence of polynomials sat- 
isfying 

t 3 u n {t)dv n {t) = 0, j = 0, ...,n-l, 
and {R n } be the sequence of corresponding functions of the second kind, i.e., 

Rn{z) := [ ^dv^t) = ^\dv n {t). (5.2) 

J z-t u n (z) J z-t 

Then, for all n large enough, the polynomials u n have exact degree n and therefore 
can be normalized to be monic. Under such a normalization it holds that 

R n n = [l + o(l)] 7 ?S n Really uniformly in (5.3) 



where S n :— S^^Kip) 71 , 7« := 2k Gy n , and k and ro were defined in (2.2 1 



Proof. This theorem is an adaptation of [Tj Thm. 3]. To see this we need several 
observations. Firstly, the orthogonality relations in [7, Thm. 3] are considered on 
Jordan arcs connecting —1 and 1, of which the interval [—1, 1] is a particular case. 
The current setting can be easily deduced by applying a linear transformation 
l(x) = [(d-c)x + d+c]/2. 

Secondly, {h n } is taken in (7} Thm. 3] to be a family of Dini-continuous 
non- vanishing functions on [c, d] such that any sequence in this family contains a 
uniformly convergent subsequence to a non-vanishing function and the moduli of 
continuity of h n are bounded by the same fixed modulus of continuity. Clearly, the 
normality of {h n } yields that all these restrictions are met in the present case. 

Thirdly, only the case m — is considered in [7, Thm. 3]. However, the 
general case we are dealing with is no different. Indeed, choose 2m zeros of each 
polynomial v n + m that converge to some fixed point in 2) (recall that the counting 
measures of zeros of v n converge in the weak* sense) and pay the polynomial, say 
P2m,ni vanishing at these points to h n . Then {h n /p2 m ,n} is again, a normal family 
of holomorphic functions with the required properties and the new polynomial 
factor v n+rn /p2 m of v n has degree no greater than 2n. 

Finally, in order to appeal to [7[ Thm. 3] , we need to show that the functions 
r„ = r2n(v n ', •) are such that r„ = o(l) locally uniformly in £), \r^ \ = 0(1) on 
[c,d], and the moduli of continuity of o ip -1 ] are bounded by the same fixed 
modulus of continuitjj^] To do so, consider 

f n {z) ^rnty-^z)) = TT T— -> E n :={^(e) : eeE n }, zeD. 
1 1 1 — cz 

Observe that the sets E n lie at fixed positive distance from [c, d] by the assumption 
supp(tf ) C Dt and therefore there exists si < 1 such that E n C D Sl for all n. Thus, 



^Observe that ip 1 is just the Joukovski transformation. 
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the Blaschkc products 



ceE„ 



z — c 



o(l) 



locally uniformly in D [H Thm. 2.2.1]. Further, 

|V>(2i) -il>(z 2 )\ < max\ip'(z)\\zi - z 2 \, z x ,z 2 G K := ^ _1 (D S J, 
and we get from the admissibility of {E n } that 

J2 \c- A„(c)| < s 2 , A„(c) := 1> (A„ {ip-\c))) , 

for all n and some positive constant s 2 . Consider now the functions 



(r n /kKz)= [] J 



1 — cz 



n r 



1 — cz 



A„(c)z' 



2 e 



Clearly, this is a sequence of outer functions in D. Moreover, 



log | (?„/&„) (r) | = lo S 



c£E n 



1 



(A„(c) - c)r 



< 



A„(c) 



< 



A - 2 



c£E„ 



A„(c)| " 



1-A„(c)r 

for r G T. Thus, we have that 

\r n \ < s 3 \b n \ = o(l), s 3 := exp{s 2 /(l - s x )} 
locally uniformly in D and 

|r„| = |r„/6 n | < s 3 on T. 

Therefore, the corresponding properties of r„ and \r^\ follow. 

Next, we show that \r^ ° ip" 1 ] have moduli of continuity majorized by the 
same function. As \b n \ = 1 on T, it is enough to consider \f n /b n \. Let tx,t 2 G T. 
Then 



(5.4) 



log 



(r n /b n )(T 2 ) 



(r n /b n )(Ti) 



1 



= E lo S 

Therefore, we have with S4 := s 2 /(l — si) 2 that 
exp{-s 4 |ri - r 2 |} < 



(n -r 2 )(c- A„(c)) 



(l-cn)(l- A„(c)r 2 ) 



< 



52 In - T 2 | 
(1- S1 )2 



(r„/6„)(n) 




(r n /b n )(T 2 ) 


(r n /6 n )(r 2 ) 




(r n /b n )(Ti) 



< exp{s 4 |n - t 2 \}. 

Moreover, denoting by Arg(z) G (— tt, it] the principal argument of z 7^ and using 
|Arg(l + z)\< arcsin \z\ < n\z\/2, \z\ < 1, 

we get that 



Arg 



(r n /b n )(r 2 ) 
(r»/b n )(n) 



Arg 1 



(n-T 2 )(c-A n (c)) 
(l-cr 1 )(l-A„(c)r 2 ) 



< 



S4 n - t 2 



2/tt 
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for s^\t\ — T2I < 1. Hence, for such n and Ti we obtain that 



(r n /b n )(Ti) - (r„/6„)(r 2 ) 



< s 3 



1 



< S 5 \Ti - T 2 \ 



(r n /b n )(n) 

for some absolute constant S5. This finishes the proof of the theorem, granted [7J 
Thm. 3]. □ 

6. Proofs of Theorems [IH4] 

We start by providing several auxiliary results. 

Lemma 8. Under the conditions of either Theorem^ or Theorem^ it holds that 
q n = u n - m q n! m, q n ,m = (1 + o(l))?, (6.1) 

locally uniformly in C \ z{q), where deg(u n - m ) — n — m and deg(q n ^ m ) = m. 
Moreover, the zeros of the polynomials q\ form an admissible interpolation scheme. 

Proof. It follows from iHJ Thm. 2.4] that if supp(/i) is a regular set with respect 
to the Dirichlet problem, ji has an argument of bounded variation on supp(/i), 
and |/i|([£E — 8,x + 5]) > 15 L for all x G [c, d) and some fixed positive constants 
I and L, then in any neighborhood of rj € z(g) the polynomials q n have at least 
m(rj) zeros for all n large enough (in fact, no more then 771(77) pl us an absolute 
constant depending only on f), which is indeed equivalent to ( |6.1[ ). Clearly, all 
these requirements on the measure fj, are met in the present case. 

Concerning the admissibility property, observe that the zeros of q n are con- 
tained in D by the very definition of g n and their counting measures converge 
weak* to the Green equilibrium distribution on [c,d] by [SJ Thm. 2.1]. Thus, the 
second requirement for admissibility is satisfied. So, it only remains to construct 
the rearrangements A„ that we shall simply take to be the identity mappings. 
This way we are required to show that the sums \4>(£,j. n ) — <ft(£,j,n)\ remain 

bounded when 77 — > 00, where £j- n are the zeros of q n and (/>(•) := tp(l/-). Since 4> 
is holomorphic in D* , {£j>} C D, and D C X>*, it holds that that 

|<^(-2i) - (f>(z 2 )\ < 8\\Z\ - z 2 |, Zi,Z 2 e D, 

where s\ := max-f |^'| = maxj |</>'| = maxn |</>'| by the very definition of cj> and the 
maximum modulus principle for analytic functions. Hence, 

n n n 

E - < 2si E |lm(0,„)l < 2*i £(7r - Angle(0,„)). (6-2) 

j'=i i=i 3=1 

where Angle(z) := |Arg(a — z) — Arg(6 — z)\, Arg(z) E (— tt,it] is the principal 
branch of the argument of z, and we set Arg(O) = n. The uniform boundedness of 



the sums on the right-hand side of (6.2) was established in [SJ Lem. 3.1], using in an 
essential manner that the argument of /i is of bounded variation, as a prerequisite 
for the proof of Thm. 2.1]. This finishes the proof of the lemma. □ 
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Proposition 9. Under the conditions of either Theorem^ or Tht 
that 



(f-ftO(z) 



0) f (b 2 n w n q)(t) dfi(t) 



Q%w n q)(z) J q n ,m(t) z-t 



z e 



lorem [#] it holds 
ns. (6.3) 



Proof. Let {«„} be a sequence of singular vectors associated to {g n } having inner- 
outer factorizations (3.51. It was obtained in [31 Prop. 9.1] that 



Kf(Vn)(Q = Vn£ (bnjnWn (£) = <J n (b n j n Wn /2 ) (£), p > 2 ^ ^ 

= ? (b n a n )(0 = (b n a n y (£), p = 2, 

for a.e. £ 6 T, where j n is some inner function and a n £ -ff 2 . Following the analysis 
in [5j Sec. 10], this leads to orthogonality relations of the form 



(l n -iq n w n )(t) 



dn(t) 



(l n -iq n Wn)(i~) p(r) dr 
SO) q(t) 2th 







(6.5) 



for any polynomial l n —\, deg(l n _i) < n — 1. In another connection, (3.3 1 yields 
that 

w n Tif(v n ) _ V-(fv n ) 



Wn(f~ 9n) 



(6.6) 



The right-hand side of ( 6.6 ) is holomorphic outside of D and is vanishing at infinity. 
So, by the Cauchy theorem it can be written as 



V-ifVr, 



<z) 



1 



K{z) 

Qn(z) 



(K)(t) dr 
z — t 2iti 
(q n Wn)(i) d(J,(t) 



1 

2tH 



{q n w n )(T) p(r) dr 



q n (t) z - t 

for \z\ > 1. Using ( 6.5 1 with Z n _i(t) = (q n (z) — q n (t))/(z — t), we get that 

-P-(fv n ) = ( r (q n w n )(t) dfxjt) + 

■ ' " ' ' ' • ' 2iti 



bn 



qn(z) 



Z - t 



Qn(r) q(r)z-T 
t), we get thai 
(q n w n )(T) p(t) dr 



q^{r) q(r)z-T 



Applying ( [675] ) again, now with Z n _i = ((qu n - m )(z) - (qu„- m )(-))/(z - •), and 
using the Cauchy integral formula to get rid of the second integral, we obtain that 



P-(fl>n) 



(Z) = 



(u n - m qqn)(z) 



qq n w n )(t) dfi(t) 



z - t 



\z\ > 1. 



Observe now that the last expression is well-defined as a meromorphic function 
everywhere in D. Thus, it follows from (6.6 1 that (6.3 1 holds. □ 

Lemma 10. Let \ n be a sequence of Borel complex measures on [c, d] such that 

q(t) d\ n (t) 



F n (z) := 



zeD, 



q n ,m(t) z-t 

converges to some function F locally uniformly in T>. Then 

d\„{t) 



F n -F = o(l), F n (z) 



zeD, 
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locally uniformly in D. 



Proof. Assume first that q(z) = (z — rf) m . Let z G Df and Li and F2 be two Jordan 
curves encompassing [c,d] and {77}, respectively, separating them from each other, 
and containing z within the unbounded components of their complements. Then 



F, 



In 



1 

2^i 



771 

r 1 \ 1 



F n {z) + 



1 




(m — 1)! \ z — ■ 

where we used the Fubini-Tonelli theorem and the Cauchy integral formula. Thus, 



F n (z) = (q n . m F n )(z) - T^ im _i(g n>m F„; z) /q(z), Z6S, 



(6.7) 



where T^ iTO _x(/; •) is the (m — 1)— st partial sum of the Taylor expansion of / 
at 77. By (6.1) the polynomials m —i(q ntm F n ; •) converge to zero as n tends to 
infinity and the claim of the lemma follows by the maximum modulus principle 
for analytic functions. By partial fraction decomposition, the case of a general q 
is no different. □ 



Proof of Theorem [7] Let 



1 Qn-\-7ji.mQV^n-\-m 7 Qn-{-'m.mQ'^n-\-m^ / '^ 7 

dv n := ^ dfi = ^ dU[ Cid ]. 



a 2 



a 2 



(6.8) 



Then we get from (6.5 1 applied with n replaced by n + m and l n +m-i(t) = t J q(t), 
j = 0, . . . , n — 1, that 



t j u n (t)dv n {t) = 0. 



(6.9) 



Hi 



So, the asymptotic behavior of u n is governed by Theorem|7] applied with v, 
and h n = q n +m.mqw n+m , due to Lemma [8] and the fact that {w n } is a normal 
family in D none of which limit points has zeros. The latter was obtained in [9, 
Lem. 3.4] under the mere assumption that \i has infinitely many points in the 
support and an argument of bounded variation. 
In another connection, observe that 



and that 
that 



b n (Wn< 72 ) " = \bn\ 2 (W /2 ) " = (j„< 72 ) " on T 
(jnWn is the trace of a function from Hq . Thus, it follows from (6.4) 

V-(b n n f (v n )) = a n V- (b n (6„j„<' /2 ) CT ) = o-n (jnwC 2 )" ■ 



It is also readily checked that 

T-(b n H f (v n ))(z) = V_(b n V.(fv n ))(z) = V-(fb n v n )(z) = 



(fb n Vn)(r) dr 
z — t lixi 
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for \z\ > 1. Hence, we derive by using the Fubini-Tonclli theorem that 



7n- 



In- 



(b 2 n w n )(t) 

z - t 



d/j,(t) 



(b 2 n w n )(T) p(r) dr 
z — t q(r) 2m 



(6.10) 

for \z\ > 1, where 7„ has the same meaning as in Theorem[7] As the right-hand side 
of (6.10) is defined everywhere in £>f, the restriction \z\ > 1 is no longer necessary. 
This, in particular, implies that j n is a finite Blaschke product as neither singular 
inner factors nor infinite Blaschke products can be extended even continuously on 
T. However, notice that first we should evaluate the second integral on the right- 
hand side of (6.101 by the residue formula and only then remove the restriction 
\z\ > 1. Clearly, this integral represents a rational function vanishing at infinity 
whose poles are those of q. It is also easy to observe that if 

l-l m {b 2 n w nP )^{7 1 ) = o{l), k = 0,...,m(r ] )-l, and (b 2 n w nP )( V ) + 0, (6.11) 

this rational function converges to zero locally uniformly in C \ z(q) and has poles 
of exact multiplicity 771(77) a ^ each r\ £ z(q). Now, we have by ( |6.3| , (5.2 1, and (5.3 1 
that 

Qn,m ^in—vaVn 1 



and 



(f-ffn) 

(t) dv n ^ m {t) 

'Jn — m % v 



blw n 



y„ 



In 



n—m J -^n—m j 



(b 2 n w n )(t) q(t) dfijt) = 1 + 0(1) 
Jn-m 1n,m(t) Z - t ro(z) 



(6.12) 



(6.13) 



locally uniformly in D. Then we get from (6.12) that 

Vn){z) = j n ^ m (blw n )(z)(f^q + p)(z) +J n l m (qb n W n fl n 
where f = f M +p/q and g n — h n /b n , and therefore for 77 € z(q) we obtain 

(q n , m y n ) {k) (v)=ln- m (b 2 n w n p)^( V ), k = 0,...,m(r 1 )-l. (6.14) 



Hence, the first part of (6.11 ) follows from ( |6.1[ ) and the normality of {y n }, which 
is immediately deduced from (6.131. The second part of (6.11 1 holds since (6.131 
and (6.14), applied with k — 0, yield that for all n large enough we have 

m \O n W n pU n —m 

and therefore (b n w n p)(rj) as well as (b^Wnp)^) cannot vanish. 



In another connection, (6.13) and Lemma 10 yield that 

(fy»n)V), l + o(l) 



1 



In- 



Z-t 



-dn{t) = 



»(*) 



(6.15) 



locally uniformly in £). Thus, combining (6.11 1 and (6.151 with (6.10), we get that 

^(jn<' 2 \ 

■n-rn. \ / 



1+0(1) , t n 



In- 



to 
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locally uniformly in 35, where deg(£ n ) < m, the polynomials £ n are coprime with 
q, and converge to zero locally uniformly in C when n — > oo. Equivalcntly we have 
that 

^-j n wg/*W = 1 + o(l) + ^ = 1 + o(l), (6.16) 

where the first o(l) holds locally uniformly in 35*, the second one holds locally 
uniformly in 35*, and £ n (z) := z m ~ 1 £„(l/z), deg(£„) < m since deg(£ n ) < m. 
This, in particular, implies that the Blaschke produ cts j n are identically 1 for all 
n large enough since the right-hand side of (6.161 cannot vanish in D for such 

n. Finally recall that by its very definition wfi^ 2 has unit L 2 norm. Therefore, 
deformation of the integral on T to [c, d] covered twice by the Cauchy integral 
formula yields that 



1 



\dr\ 
2ir J T |tt)(r)| 



o(l) 



1 



dr 



1 



dt 

cA (n»-fi)(t) 



2ni J T (roro)(r) 
o(l)=T- 2 + (l) 



o(l) 



(6.17) 



since t o(r) = rro(r) on T, K> = — tn + = — i\tt> | on [c, d], and on account of (2.151. 
Thus, (3.6) shall follow from (6.161 and (6.171 with l n := T£ n upon showing that 

7»/l7«l = l + o(l). (6.18) 
The latter is an easy consequence of (6.161 since tr(0) = 1 and w n (0) > 0. □ 



In the next proposition of technical nature, we define a special sequence of 
Szego functions for the condenser 35 fl 35* that appears in Theorem [5J 

Proposition 11. For each p € [2, oo] there exists a normal family of non-vanishing 
functions in 35 n 35*, denoted by {D n }, such that 



M\Dl 



{r n /r)±S± 



Wf.1- 



z,d), 



(rD n ){l) 



>0, 



where r n := r n (q n ; ■), r := r m (q; ■), M 2 = M/\G b 2 wfl \, and M is given by ( |3. 10 > . 
Moreover, each D n satisfies D n (z)D n (\/z) = 1, z € 35H35*, has continuous traces 
on each side of [c,d], and has winding number zero on any curve separating [c,d] 
from [c, d]* . 

Proof. The concept of Szego function for a condenser initially was developed in [TSJ 
Thm. 1.6] in the case of an annulus. It was shown that if Y is a continuous (strictly) 
positive function on T s , s < 1, then there exists a function u, harmonic in A s , such 
that u — |log(y/My) on T s , u = on T, and u(l/z) = —u(z), z 6 A s , where 

My '■= exp | j T log Y(t) j is the geometric mean of Y. Moreover, it was shown 

that u has single-valued harmonic conjugate v. Moreover, the latter is unique up 
to an additive constant. Finally, it was deduced that D S (Y; •) := exp{w + iv}, the 
Szego function of Y for A s , is a non- vanishing holomorphic function in A s such 



20 



M. Yattselev 



that M Y \D s (Y;-)\ 2 = Y on T s , D S (Y; z)D s (Y;l/z) = 1, z € A s , and D S (Y;-) is 
an outer function in A s with zero winding number on any curve in A s . The latter 
was not explicitly stated in [TB] but clearly holds since log \D S (Y; -)| = u and 
therefore it is the integral of its boundary values against the harmonic measure 
on dA s while Arg(D s (Y; •)) = v, which has zero increment on any curve in A s . 
Obviously, the Szego function for A s is unique up to a multiplicative unimodular 
constant. 

Let now y + and y~ be two continuous positive functions on [c, d] whose values 
at the endpoints coincide. We define the geometric mean and the Szego function 
of the pair rj := (y + ,y~) for the condenser 3nS* by 

M„ := My and D v (z) := D p (Y; (p{z)), z e £ n £>*, 

respectively, where Y(t) — y ± (f^ 1 (r)) , r S . It is an immediate consequence of 
the corresponding properties of D P (Y; •) and tp that is outer, has non-tangential 
continuous boundary values on both sides of [c, d] and [c, d]* whenever X) is a Dim- 
continuous paiiQ has winding number zero on any curve in S (~l S)*, and satisfies 



M„/y±(lA), te[c,d]*, 



and 

Now, put y : 



l, z e D ns*. 

|& 2 u>/i| and 



j/+(t)j,-(t) = IS+^WS^ 



(*)| : 



(& 2 u;/i)(i) 



G b 2 t 



Observe that in this case 
\ y+(t)y-(t) ^ 

\Gb 2 wu\ 2 



t £ [c,d], by (2.7) and (2.9 1. Then for t)„ := {yt,Vn) we § et tnat 



logM„„ = 



io g [r w (r)y n (f)]^i 



log[y+(t)y"W] 



|dr| 



log 



y + (t)y-(t) 

\Gb 2 wfi\ 2 



\*[\ 
2wp 



2np 
logMt, — log 



where £ = ip 1 (r) = 95 (t). Further, by (2.151 and (3.10), we have that 



log M„ = 



c.d] 



7T/9 



So, Mf, n = A//|Gf,2 u , /i | = M 2 and therefore the claim of the proposition follows by 
setting D n := D^ n with the chosen normalization (recall that the functions D^ n 
are uniquely defined up to a unimodular constant). □ 



The following Lemma was proved in [SI Lem. 4.7]. 



4 This means that Y and therefore logy are Dini-continuous as tp 1 is Lipschitz on T p . Hence, 
the boundary values of v axe continuous on 8A P 1121 Ch. III]. 
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Lemma 12. Let U be a domain, dU = K1UK2, K\ and K2 be two disjoint compact 
sets in C, and u be a harmonic function in U . Assume that 

du 



On 



ds = 0, 



where d/dn and ds are, respectively, the normal derivative and the arclength dif- 
ferential on V , and the latter is an oriented smooth Jordan curve that separates 
K\ from K2, has winding number 1 with respect to any point of K±, and winding 
number with respect to any point of K2 ■ Then 

sup limsup u{z) > inf liminf u(z) 
z'&Ki z->z' , zEU z'£K 2 z->z', zeu 

and the same relation holds with K± and K2 interchanged. 

Proof of Theorem^ It was shown in Lemma[8]that q n can be written 

and the behavior of u n - m is governed by Theorem [7] with v n defined in (6.8 ). Thus, 

we have from ( |5.3| that 

1 + o(l) = u n - m S n - m = u n - m (Ki()) n ~ m S q p,/g2S qn mWn (6.19) 

locally uniformly in D. Further, since q n ,mW n converges to qw uniformly on [c,d\ 
by Lemma[8]and Theoremjljwe get that S qn mWn = [1 + o(l)]S qw uniformly in £) 
and therefore we obtain from ( |6.19 ) that 



'S, 



1 + o(l) = 1t„_ m K™ ' 

locally uniformly in D, where b 2 — q 2 jq 2 . Now, it follows from (2.101 that 

ip m - n q r n 



(6.20) 



S. 



'i 2 hi 



_ q2 _ 



Gq/q„ Qn 



where r := r m (q; •) and r n :— r n (q n ;-) as in Proposition 11 Hence, we deduce 
from ( |6.20[ ) that 

b n q K n - ,n r^ 



0(1) 



b q 



locally uniformly in D. Then Lemma [8] implies that 

\ n X n D n b n v m - n b- 1 = l + o(l), 
locally uniformly in S)t, where 

K ■= — „ and X n [z) := — — - — 

Gghr, MD n r 

Now, we shall show that 

A n ^n = l + o(l) uniformly in B 

Observe, that 



(6.21) 
(6.22) 



\X±\ = 



( Sv 2 



w/j, 1 n 



\MD n 



1 on [c, d] 



(6.23) 
(6.24) 
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by the very definition of D n (see Proposition 111. Moreover, since the zeros of 
r, z(§), lie outside of D and the zeros of r„, z(g„), approach [c, d]* and z(g) by 
Lemma [8] and Theorem [7J the functions X n are zero free in some neighborhood of 
D, where the values on [c, d] are twofold. Further, the winding number of X n along 
any smooth Jordan curve encompassing [c, d] in D is equal to zero. Indeed, the 
winding number of Sf,2 w f l /D n on such a curve is zero by the properties of Szego 
functions, r n /r has winding number m — n since it is meromorphic outside of [c, d] 
with n zeros and m poles outside of D, and it follows from [18j Ch. VI] that ip has 
winding number one on any such curve. Thus, logX n are well-defined holomorphic 
functions in D \ [c, d]. In turn, this means that log \X n \ satisfies the conditions of 
Lemma 12 with U = D\ [c, d). Applying this lemma in both directions, we get 
from (16.24} that 



inf|X n |<sup|X±| = l = inf \X±\ 
T [c,d] Ml 



< sup \X n \ 

T 



In another connection, (|6.21} and (|6.1} yield that uniformly on T we have 

o(l) 



\X„ 



G q/q n 



M(np) n - m Gq/q n 



6^ W LL 



1 



(6.25) 



(6.26) 



since D n , ip, b n , and b are unimodular on T. Combining (6.25) with (6.261, we get 
that |A„| = 1 



o(l) and therefore 

(nip) n ~ m S b 2 w(t r n 



|A„X„|= ^ — — =l + o(l) uniformly in D (6.27) 

by the maximum principle for harmonic functions applied to ± log \X n | in D\ [c, d] . 
Hence, {X n X n } is a normal family of harmonic functions in D \ [c, d] and all the 
limit points of this family are the unimodular constants. Therefore (6.231 follows 
from the normalization of D n (see Proposition 11 1 and the fact that <p(l)=l. 
Clearly, we can rewrite ( 6. 22} with the help of (6.231 as 



D n b n <p r 



n b 1 = 1 + o(l) uniformly on compact subsets of DnBf 



Now, recall that D n (l/-) — 1/D n . Moreover, the same property holds for b n , 
and (p. Thus, 

D n b n ^ m - n b- x = 
uniformly on closed subsets of 23? \ ! 



-»6" 1 )(1A) = l + o(l) 



and (|3.7} follows. 

By the very definition of j v 



\ U q/q n 



It only remains to prove (3.8) and ( |3.9 
Theorem |7] we have that 

IT"-™! = K 2(n-m) = t 1 + °( 1 )J ^-m) 

by Lemma |8] and limit (3.6 1. Further, the very definitions of M and A„ yield that 



2\Gb2 l 



2\G, 



\M 2 



\G 



q/qn 



2Mp 



2{n-m) 



^2(n— m) 



A,. 
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Since |A n | = 1 + o(l), it holds that 

7n _ m = [1 + o(l)]|7„-m| = [1 + o(l)]2Mp 2( - n - m \ (6.28) 

where we used (6.18). Thus, ( |3.8| follows from ( |6.17 1. Finally, we deduce from 
( |6.12| and ( |6.13| ) that 

(f-<7„) = [l + o(l)] 



1 + 0(1) 



q 62 w„ ro 

uniformly on compact subsets of IDf PI D. Since q n ,m/l 1 by Lemma|8] ic„ — * w 
by ( pL6) , and using ( |6.28| ) with ([3J|, (|3J]) follows. □ 

Proof of Theorem^ Recall that Lemma [8] holds under the conditions of this the- 
orem and the Szego functions D n exist for p = 2 as well. As the starting point of 
the proof of Theorem [2] was the application of ( |5.3| with v n defined in (6.8 1, all we 
need to do is to show that the conditions of Theorem [7] still hold under the present 
assumptions. This is tantamount to show that T x , defined in (5.1l, is minorized 
by defined in the statement of the theorem, for all x € x. In other words, that 



a; € x. 



ty x < lim inf min { I (r n Q„ 
Equivalently, we need to show that 

— liminf log |r,^(a;)| 2 — lim sup log |rj(x), 

n— >oo " n->oo " 1 — So 

by Lemma [8] and the definition of Let, as in the proof of Lemma [8] £j 
j = 1, . . . , n, be the zeros of q n . Then we get from (6.2 1 that 



2 4s 1 (y h +2m7r) 



x e x, 



X>(V&*)-^1/&>)|<2*1^ 



(tt - Anglefe- „)) < 2si(14 + 2m7r), 

:= r„ o Then, 



where we used [51 Lem. 3.2] for the last inequality. Put 
exactly as we did to prove (5.4), we obtain that 

log |f„| < 2si(T4 + 2m7r)/(l - «„), 

which finishes the proof of the theorem. □ 

To prove Theorem [4] we need the following lemma. 

Lemma 13. Let R = P/Q be a rational function of degree d, ( G C. and S > 0. 
Assume further that P and Q have no zeros in {z : \z — £| < 8}. Then for any 
k < d, k GN, there exists Ck = Ck{8) independent of R such that 

R (k) (C)/R(C)\<c k d k . 

Proof. Clearly, if T is a polynomial of degree at most d with no zeros in {z : 
\z - CI < 5}, then 



tO')(c) 

no 



< 



d-...-(d-j+i) 

8i 



< 



d> 3 
5. 



J = 1, 
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Thus, it can be checked that 

3 



1 \ W 



tW(0 



di 



<C*d\ 



j = 1, . . . , k, where coefficients cj,{di},{*<} d° not depend on T. Then 



i?( fc )(C) 



QiO 



Q(0 



(k-j) 



□ 



Proof of Theorem^ As the forthcoming analysis is local around 77, we may sup- 
pose without loss of generality that m(rj) — m, i.e. r\ is the only zero of q. 
Exactly as in ( |6.14 ), we obtain that 

2 

y^(v) = (Y n q n ,m) {k \v), fc = 0,...,m-l, Y n := (6.29) 

'Yn—mQn 



It is apparent from (6.221, (6.231, and foremost (6.211, which holds locally uni- 



formly in 35 rather then S)f, that 

— rn. 



1 + o(l) locally uniformly in DflS. 



So, we see using (|3.6|), <\6.1\, (|6.17|) with (|6.18), and (|3.8|) that 

•o(l) 



F n = [1 + o(l)] 



^2(n-m) 



1 

2M 



2(n— m) 
if \ wp 



Dlq 2 



(6.30) 



uniformly in some neighborhood of rj. Thus, we obtain from (6.29) with k = that 

Xn1n, m (v) = 'I, Xn ■= ~Y n (v) / Vn(v) , (6-31) 

where for each n we fixed an arbitrary root Xn- Observe also that Xn tends to infin- 



ity geometrically fast by (6.301 since > p and the boundedness of {|y n (??)|}, 



which is apparent from (6.13). By putting k = 1 in (6.291, we see that 

i (y 7 'M y'nivY 



Y m "Y 
An yn.m 



Xn 



0(1) 



since {v^iif) / y n (v)} 1S a convergent sequence by (6.131 and Y n are rational func- 
tions, which do not vanish in some fixed neighborhood of r/, multiplied by w n , which 



form a convergent sequence by (3.6l, the numbers K^/^nX 7 ?)! grow linearly with 



n by Lemma 13 while 1/xn decays exponentially. Continuing by induction, we get 

o(l), (6.32) 



v m-.( fc ) M ^ fk\ Y^( V ) X%- k+j q ( n,™ j) (v) 1 2/1%) 

Xn Qn,m\ I ) 



3- 



Xu 



Xn Vn(v) 
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for any k = 2, . . . , m — 1. Hence, we deduce from (6.31 1 and (6.32 1 that 

m m— 1 

n (z + x»fa - %,„)) =z m +J2 Xn- k qi%(v)z k = z m + 0(1) - 1, 

k=l k=0 

uniformly in some neighborhood of rj. In particular, this means that 



, l+<5fc,n 

.„ // • exp 



f 27r/ci 



V 



1. 



where Sk, n = o(l) for each fc and IlfcLi(l + ^fe.n) = 1- By setting 



/ \ ' 



i + s, 



2(m — n) 



= -[l + *fc,n] 



2(m— n) 



we see that (3.11 ) follows. The boundedness of {max^ \ AZ |} is a consequence of 
d6~30l) and doST □ 



7. Proofs of Theorems |5] and H] 

Proof of Theorem [5[ Let q„ be the denominators of II„ . We start by showing that 

q n = u n - m q ntm , q n , m = (1 + o(l))q, (7.1) 

locally uniformly in C \ z(q). This follows from jTUJ Thm. 2.4] in the same fashion 
as (6.1 1 followed from [5J Thm. 2.4]. The requirements placed on fi are the same, 
so they are satisfied. However, in [TUJ Thm. 2.4] there are also restrictions placed 
on the interpolation schemes. Namely, an interpolation scheme S should be such 
that supp(<?) n ([c, d] U z(q)) = 0, the probability counting measures of points in 
E n would converge to some Borel measure with finite logarithmic energy, and the 
argument functions of polynomials v n would have uniformly bounded derivatives 
on [c, d] . 

Clearly, the first two requirement placed on the interpolation scheme is the 
second requirement of the admissibility property. Hence, we only need to show the 
uniform boundedness of the derivatives of the arguments of v n . Clearly, it amounts 
to show that 



lim sup ; 



1 



Im 



lim sup ; 



1 



:,d\ 



E 



Im 



1 



< oo. 



(7.2) 



led] 



Since Im(f — e) = Im(e 
1 



t) for t £ [c, d] , we have that 



Im 



E Im 

E 



Im 



A n (e) 
A n (e) 



E Im 



(t - A n (e))(t - e) 



< 



-y 

2s 2 



|A„(e) 



where the sums are taken over e £ E n n C and s > is such that \t — e\ > s for all 
e £ E n and n £ N. So, (7.2 1 and therefore (7.1 ) follow from the admissibility of S . 
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It is well-known |25l Lem. 6.1.2] and is easily seen from the denning properties 
of Pade approximants and the Fubini-Tonclli theorem that 



t 3 q(t)q n (t) 



dfx(t) 



and 



V n (t) 

v n (z) 



0, j = 0, ...,n- m- 1, 

n — rn )(*) Mt) 



(q n ql n - m )(z) 



z - t 



V n (t) ' 



(7.3) 

z G % (7.4) 



for any polynomial „_ m Z of degree at most n — m. Now, using decomposition ( |7.1| 
and denoting 



:= 1 — a/i = 1 awr, 



orthogonality relations ( |7.3[ ) become 

t j u n (t)dv n (t) = Q, j = 0, . 



(7.5) 



, n — 1. 



It is also quite easy to see that the asymptotic behavior of u n is governed by 
Theorem [7] applied with /i„ = q n + m . m q. The orthogonality relation above also 
imply that 



R n (z) := 



U n (t) 1 

-^„(i) = — ^ 



t 



dv n (t), z eD. 



Thus, putting l n - m = i„_ m , we can rewrite (7.4 1 as 

I = • 

Hn—mqn.mq 



Hence, we derive from (5.3) and (7.1l that 



f -n n = [i + o(i)] - 



n>qq n .r 



[2 + o(l)]G i , n _ m Si n 



K>q n . m q 



(7.6) 



(7.7) 



(7.8) 



locally uniformly in D. Therefore, we get from (|2.10| and (|7.5|) that 



qn,m,q 



r 



G q S VnTj) 2 



q n , m ?P m q^ m G Vn Sl n 



(QVi,m ! ' )^ 



[l + o(l)]G fl S i 



2 ' n 



(7.9) 



locally uniformly in 2)j, where r n and r are defined as in the statement of this 
theorem. Combining (7.8 1 with (7.9 1 we get (4.2 1. 

Finally, observe that the boundedness of the variation of argument of h was 
needed in order to appeal to |10l Thm. 2.4]. However, when the rational summand 
of f is not present (q = 1), (4.2 1 is a consequence of Theorem [7] only and the latter 
does not require the boundedness of the variation of argument of h. □ 
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Proof of Theorem^ The asymptotic equality in (4.3 1 is exactly the one in (7.1|. 
The fact that u n - m have no zeros on compact sets in D follows since the asymptotic 
behavior of u n _ m is governed by Theorem [7] with v n given by (7.5 1 and all the 
zeros of such orthogonal polynomials approach [c, d] . 

Let rj € z(g). As in the proof of Theorem [4] we may suppose without loss 
of generality that m(rj) = m. i.e. // is the only zero of q. Using the notation of 
Theorem [7] we can rewrite (7.7 1 as 



or equivalcntly 



Vn 



lnV n 



(f II n ) U n — rn R n - 
^Jn—rn 



in (UQn - Pn) 

"^in—mVn 



Y r , := 



7 



(7.10) 



n—m u n 



where f = f M + p/q. It follows from (5.3 1 that 

y n W = 1 + o(l) locally uniformly in T>. 
In particular, it means that sequences {|yn^(f?)|} are uniformly bounded above 



and away from zero for all k € N. Moreover, (7.10) yields that 
y { nHv) = (Y n q n ^ k H V ), k = 0, . .. ,m - 1. 
This, for instance, implies that neither of r]k,n, k — 1, . . . , m, the zeros of q n , m , is 



equal to 77. Further, using (5.3l, (2.10), and (7.1 1, we get that 



Y n = [1/2 + (l)]V 2 >/(G Ag2 5? 92 r- n ) 

uniformly in some neighborhood of rj. Now, it it clear that we may proceed exactly 
as in the proof of Theorem [3] with the only difference that we set 
/ \ m 

:= ( X + S k,n)r n (r))/X™ - -(1 + 6 k , n )(r n y n /Y n )(r,). 

□ 



8. Numerical Experiments 

The Hankel operator Ttf with symbol / e H°° + C(T) is of finite rank if and only 
if / is a rational function [THl Thm. 3.11]. In practice one can only compute with 
finite rank operators, due to the necessity of ordering the singular values, so a 
preliminary rational approximation to / is needed when the latter is not rational. 
One way to handle this problem is to truncate the Fourier series of / at some 
high order N. This provides us with a rational function /jy that approximates / in 
the Wiener norm which, in particular, dominates any L p norm on the unit circle, 
p € [1 , 00] . It was proved in |15j that the best approximation operator from 



(mapping / to g n according to ( 3 .3 1 ) is continuous in the Wiener norm provided 
(n+l)-st singular value of the Hankel operator is simple. It was shown in [21 Cor. 2] 
that the last assertion is satisfied for Hankel operators with symbols in some open 
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dense subset of H °° + C (T) , and the same technique can be used to prove that it 



is also the case for the particular subclass (1.1 1. Thus, even though the simplicity 
of singular values cannot be asserted beforehand, it is generically true. When it 
prevails, one can approximates /at instead of / and get a close approximation to g n 
when N is large enough. This amounts to perform the singular value decomposition 
of H fN (see [Ml Ch. 16]). 

As to Pade approximants, we restricted ourselves to the classical case and we 



constructed their denominators by solving the orthogonality relations (7.3 1 with 
v n = 1. Thus, finding these denominators amounts to solving a system of linear 
equations whose coefficients are obtained from the moments of the measure //. 

The following computations were carried with MAPLE 8 software using 35 
digits precision. On the figures the solid line stands for the support of the measure 
and circles denote the poles of the correspondent approximants. The approximated 
function is given by the formula 



m = 



7e lt dt r it+1 dt 

[-0.7,0] * - t y/(t + 0.7)(0A-t) i[o,o.4] z - t y/(t + 0.7)(0A-t) 



5!(z-0.7-0.2i) £ 




FIGURE 1. Poles of Pade (left) and AAK (right) approximants of 
degree 10. 



Acknowledgment. / express my sincere gratitude to Dr. L. Baratchart for valuable 
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FIGURE 2. Poles of Pade (left) and AAK (right) approximants of 
degree 20. 




0.701 0.699 



FIGURE 3. Poles of Pade (left) and AAK (right) approximants of 
degrees 21-33 lying in an neighborhood of the polar singularity. 
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